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^ : Abstract 

We present a new method for the statistical investigation into the distributions of the 
angle /3 between the radio axis and the normal to the galactic disk for a sample of Seyfert 
^ . galaxies. We discuss how further observations of the sample galaxies can strengthen the 

^ \ conclusions. Our data are consistent with the hypothesis that AGN jets are oriented ran- 

domly in space, independent of the position of the plane of the galaxy. By making the 
G^ ■ simple assumption that the Standard Model of AGN holds, with a universal opening angle 

^ ■ of the thick torus of (/)c, we demonstrate a statistical method to obtain an estimate of 0c- 

Q> . Our data are not consistent with the simple-minded idea that Seyfert Is and Seyfert 2s are 

^ \ differentiated solely by whether or not our line of sight lies within some fixed angle of the 

Q-i| jet axis. Our result is significant on the 2a level and can thus be considered only suggestive, 

P ■ not conclusive. A complete sample of Seyfert galaxies selected on an isotropic property is 

required to obtain a conclusive result. 



1. Introduction 



^ \ We concern ourselves here with the angle /3 at which the central accretion disk of an 

\-< \ active galactic nucleus (as evidenced by the direction of the central radio jet) lies relative to 

the host-galaxy disk, and with the angle which is the angle between the line of sight and 
the radio axis (see Figure 1). 

In the standard picture of active galactic nuclei (Antonucci 1993), it is the angle which 
determines whether a galaxy appears to us as a Seyfert 1 or as a Seyfert 2. The central 
nucleus is surrounded by a physically thick circular torus of molecular gas which ensures 
that the very central regions of the accretion flow (the broad line region) are visible only 
when the line of sight and the central disk axis (assumed to align with the radio structure) lie 
within about 0c ~ 30° of each other (Osterbrock and Shaw 1988). It is a central prediction 
of the model, in its simplest form, that values of for Seyfert Is should be on average smaller 
than values of for Seyfert 2s. The estimate of 0c ~ 30° given by Osterbrock and Shaw 
comes from an estimate of the relative numbers of Seyfert 1 and Seyfert 2 galaxies per unit 
volume of space. 
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The angle (3 serves as a link with the innermost workings of the active galaxy - either 
as the memory of the activity-provoking event, or as a diagnostic for the structure of the 
accretion disk. We develop a method for evaluating the most likely distribution of such 
angles based on an observed set of position angles as seen on the sky plane, and apply our 
method to a data set culled from the literature. This study is a precursor to a larger study 
based on complete, unbiased samples of spiral active galaxies selected on nearly isotropic 
properties. 

The minimum components for activity in the nuclei of galaxies are a central black hole 
plus a fueling source which feeds material onto the black hole (BH). The processes for 
bringing material close to the core (within the innermost lOpc) of a galaxy, either for the 
initial formation of the BH or to provide fuel to the BH, are of two sorts - those that trigger 
the gas internal to the galaxy to collapse to the core, and those that trigger gas external to 
the galaxy to collapse to the core. Internal triggers include cooling flows, bars in galaxies, 
which funnel material to the core (Shlosman et al. 1989, Thronson et al. 1989), or fly- 
bys of companion galaxies (Byrd et al. 1986) which shock the ISM of the host galaxy into 
collapsing (perhaps aided by the formation of a bar) . External triggers are mergers between 
galaxies which bring material to the core from outside the galaxy (Hernquist 1989, Disney 
et al. 1995). 

Spiral galaxies have ample cold gas available to be shocked into collapse to the core 
by fly-bys or bars, so that internal triggers are viable possibilities for causing activity in 
spirals. But spirals have so much gas in the disk that the gas can mask or interfere with 
the evidence for a recent merger. Merger models show that the gas merging with the spiral 
galaxy at random angles dissipates the angular momentum that is not aligned with the galaxy 
angular momentum, and settles into the plane of the disk so as to be indistinguishable from 
the native gas (see Walker, Mihos, & Hernquist 1996, and Mihos & Hernquist 1995 and 
references therein). 

AGN in spiral galaxies appear to show precious little correlation between the direction 
of inner jets/ionization cones and the spiral disk axis (Ulvestad & Wilson 1984; Brindle et 
al. 1990; Baum, O'Dea, de Bruyn, & Pedlar 1993; Schmitt et al. 1997). They also show 
a tendency for the broad line Seyfert Is to point out of the galaxy disk while the narrow 
line Seyfert 2s point closer to the plane of the disk (Schmitt & Kinney 1996). This finding 
may, however, not be due to the intrinsic distribution of relative position angle, but rather 
due to the manner in which the active galaxies are detected for inclusion into samples. Of 
particular interest is the zone of avoidance discovered by Schmitt, Kinney, Storchi-Bergmann, 
& Antonucci 1997 who concluded that there is apparently a zone within a cone of half angle 
20" of the spiral disk axis within which no radio axis falls. The active nucleus apparently 
knows where the galaxy plane is and avoids close alignment with it. Given that fueling by 
both internal gas and by external gas favors co-alignment between AGN nucleus and galaxy 
disk, the observed mis-alignment is intriguing and requires a theoretical explanation, it 
might indicate, for example that recently suggested ideas about radiation induced warping 
of accretion disks (Pringle 1996, 1997; Maloney, Begelman & Pringle 1996) come into play 
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both to randomize the directionahty and to produce the ionization cones. 

In this paper wc undertake a more direct statistical analysis than has been attempted 
hitherto in order to try to determine the distributions of the angle f3 between the jet axis and 
the normal to the galaxy disk. The main aim of the paper is to demonstrate the feasibility of 
such an analysis, and to elucidate what further observations might be required (apart from 
simply acquiring a larger and more uniform sample) in order to strengthen the conclusions. 
In Section 2 we describe the data set that we use. In section 3 we set up a description of the 
geometry of the situation and present our statistical method for obtaining the f3 distribution 
together with some preliminary results. In Section 4 we present the conclusions of the 
analysis and discuss both how further observations could be used to enhance the results and 
the possible implications for theoretical models. 

2. The Data Set 

We base our analysis on the data set collected by Schmitt et al. (1997, henceforth SKSA). 
The full list of objects is presented in Table 1 of that paper. The sample includes all Seyfert 
galaxies in the literature which have been resolved in high resolution radio maps and which 
show linear radio structure as defined by Ulvestad and Wilson (1984). The position angles 
of the extended radio structure were taken from the literature, while the position angles 
of the host galaxy major axis were all verified using the Digitized Sky Survey Plates and 
assuming that the galaxy disk is intrinsically circular. No galaxies in the sample lie within 
10 degrees of edge on. Such objects are apparently difficult to detect as Seyfert galaxies, 
perhaps because in such a case the active component would be entirely absorbed by the 
host galaxy disk. Our method requires a measurement of both i, the galaxy inclination, and 
the radio position angle. Since 4 of the objects in the Schmitt et al. sample of 46 have no 
measurement for i, we include 42 objects in our analysis. Of these 14 are Seyfert Is and 28 
are Seyfert 2s. 

When using a sample of objects assembled from the literature, one must always consider 
selection effects. SKSA show that the Seyfert Is and Seyfert 2s in the sample have comparable 
6cm radio fiuxes. Likewise, the Seyfert Is and the Seyfert 2s are hosted by galaxies of 
similar morphological type. However, Ulvestad and Wilson (1984) note that overall, Seyfert 
1 galaxies which are resolved in the radio display shorter elongated radio structures than 
radio resolved Seyfert 2 galaxies. The standard model of AGN (Antonucci 1993) would 
interpret this as being due to the Seyfert 1 galaxies being those which observed within 
about 30 degrees of the torus axis in order that the broad line region can be observed. The 
face-on orientation would therefore tend to foreshorten the radio emission on the sky. But 
meanwhile, such foreshortening may cause Seyfert 1 galaxies to fall out of the sample, since 
the radio emission must be linear to be included in the sample. This is a sort of blind spot 
in the sample, but it should be a randomly distributed blind spot since the Seyfert Is that 
are excluded can still have their radio jet extending from any angle out of the host galaxy 
(except for in the plane of the host galaxy - but those objects do not seem to make it into 
the sample in great numbers). The exclusion of Seyfert Is with exactly face-on tori from 
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the sample should affect absolute numbers of Seyfert I's versus Seyfert 2s, but it should not 
affect the distributions of position angle overall. Thus this might be expected to affect the 
(p distributions for Seyfert Is, but should have little impact on the j3 distribution. 

Radio flux and size generally correlate with far infrared power (Sanders et al. 1989). 
This sample contains the objects among the Seyfert galaxies which are detected and resolved 
in the radio. Thus these are also the brightest and most extended in the radio. The main 
selection effect is therefore the relatively benign one of selecting the top end of the luminosity 
function. The phenomenology of the most luminous of the Seyfert galaxies may differ from 
the phenomonlogy of the less luminous. An investigation of such luminosity effects could 
well be of interest, but will have to wait for deeper and higher resolution radio data. 

We note that, although we are using the same sample as SKSA. our approach to the data 
is different. SKSA, in common with previous authors (for example, Ulvestad & Wilson 1984, 
Baum et al. 1993) made use of just the distribution of position angle 5 between the radio 
axis and the apparent galaxy major axis. Here we make use of the fact that for almost all 
of the sample we know the galaxy inclination (assuming that the galaxy disk is intrinsically 
circular) as well as 5. This enables us to undertake a more direct approach by constraining 
the geometry of individual objects, as we describe below. To take this point further we 
show in Figure 2 the distributions of pairs (5, i) for both the Seyfert Is and the Seyfert 2s. 
Evident from the Figure the selection effect on inclination i, in that near face-on and near 
edge-on galaxies have been excluded form the sample. Using a KS test, the distributions in 
i are not found to differ between the two samples, and indeed do not differ significantly from 
a uniform distribution over the sphere within the apparent sample limits. However the 5 
distributions of the Seyfert Is and the Seyfert 2s do differ, with the probability that the two 
samples are drawn from the same distribution being less that 1 per cent. This difference can 
be seen clearly in the Figure, with the Seyfert 2s predominating in the range 0° < 5 < 30°. 

3. Statistical Analysis 
3.1 Geometry 

For each galaxy we can determine two observational parameters, i and 5. The angle i is 
the inclination of the plane of the galaxy to the plane of the sky, or equivalently the angle 
between the line of sight and the vector normal to the galaxy plane. The angle i hes in the 
range 0° < i < 90°. We use a Cartesian coordinate system OXYZ (see Fig. 1, cf. Fig. 5 
of SKSA) so that OX lies along the apparent major axis of the galaxy disk, OY lies along 
the apparent minor axis, and thus OZ is the vector normal to the galaxy plane. In these 
coodinates the unit vector in the direction of the line of sight is 

ks = (0, — sini, cosi). (3.1.1) 

The angle 5 corresponds to the position angle between the apparent major axis of the 
galaxy and the direction of the radio jet projected onto the plane of the sky. By convention 
5 is taken to he in the range Q° < 5 < 90°. 
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For a given value of 5 and i, the direction of the jet, which we denote by a unit vector 
kj is determined to he on a great circle drawn on a unit sphere centred at the origin of our 
coordinate system. In the OXYZ coordinates described above the great circle is the set of 
points: 

~ i^jXj kjy, kjz) 

— (cos S sin (f), 
sin 5 cos i sin — sin i cos 0, 

sin 5 sin i sin + cos i cos 0) , (3.1.2) 

where (j) is the angle between the vectors and k^ and formally lies in the range —180" < 
cj) < 180°. We should also note that there is a mirror symmetry to the problem about the 
apparent minor axis of the galaxy, that is about the OYZ plane. In terms of our coordinates 
this translates into the statement that reversing the direction of the OX axis leaves the 
problem unchanged. Thus, formally, the sign of kjx is not an observationally meaningful 
quantity, or in other words the jet vector in fact lies on one of two great circles which are 
reflections of each other in the OYZ plane. We have therefore simplified the discussion by 
considering just one of these great circles. 

Of course not every point on the great circle is physically accessible for the reason that 
we assume that the observable jet is the one that lies above the disk plane (as seen from 
Earth). If we define /3 as the angle the jet vector k^ makes with the disc normal OZ, then 
we see from Equation (3.1.2) that 

cos P = kjz 

= sin f) sin i sin cf) + cos i cos 0. (3.1.3) 

Thus the only relevant values of are those which give 0" < f3 < 90°, or cos /5 > 0. In terms 
of 0, this means that (f) lies in the range 0i < < 0i + 180°, where 0i is the value of 

01 = tan~^(— coti/sin5), (3.1.4) 

which hes in the range —90° < 0i < 0°. We note that physically, if < 0, then the jet 

vector is projected against the half of the galaxy disk which is nearest to us, whereas > 
corresponds to the jet being projected against the half of the galaxy disk which is furthest 
from us. From the current data we cannot discriminate between these two possibilities, but 
see, however, the discussion in Section 5. 

3.2 The minimum value of P 
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We have seen that a given pair of values of i and 5 constrains the jet vector to he on a 
great circle kj(</)) on the unit sphere centred on the galaxy. The great circle passes through 
the line of sight and makes an angle b with the plane of the galaxy (the OXY plane), where 

6 = cos~^(cos5sini). (3.2.1) 

This means that we can set a lower limit to the value of (3, which corresponds to the place 
where the great circle comes closest to the pole OZ. This value, which we denote by is 
given by 

Pmin = cos-^{sm^6sm^i + cos^iy/^. (3.2.2) 

We note, however, that if we were able to distinguish between the near and far sides of 
the galaxy, that is between +0 and —(f). then wc would be able to decide which of the two 
segments of the great circle, on either side oi (p = 0, the jet direction corresponds to. In that 
case we would have 

Pmin — cos"^(sin^ 5 siu^ i + cos^ iy^"^, 
for > 

= i,for0<O. (3.2.3) 

Furthermore, we also note that if we were able to determine if the jet is approaching us or 
receding from us (that is we were able to determine the sign of k^.k^ = cos0), then this 
would further restrict the allowable range of (p, (in particular to decide whether (f) is greater 
or less than 90°) and so enable a more accurate determination of /S^m- 

However, for the purposes of this paper, in which the observational data does not tell us 
which side of the galaxy disk is closer to us, and in which we do not know whether the jet 
is approaching or receding, wc must take f3„nn as given by Equation (3.2.2). In Figure 3 we 
display the values of Pmin for the galaxies in our sample, distinguishing between the Seyfert 
Is and the Seyfert 2s. 

We stress that while these results only give us one sided limits to /3, the results themselves 
depend only on the geometry of the situation, and do not require any statistical modelling 
or additional assumptions. In this sense the conclusions drawn from them are robust and 
irrefutable. Thus we can state that, for example, more than 80 percent of the Sy2s in our 
sample have /3 greater than 20°, more than 50 percent of them have P greater than 30°, and 
more than 20 percent of them have f3 more than 54°. For the Syls in our sample we can 
state that more than 80 per cent of them have j3 greater than 17°, more than 50 percent 
of them have f3 greater than 25°, and more than 20 per cent have j3 more than 32°. These 
findings agree with the conclusions of Schmitt et a/. (1997) who present evidence for a 20° 
zone of avoidance around the galactic pole. 
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It might be tempting to draw the conclusion from Fig. 2 that the distributions for Pmin 
(or equivalently the distributions for the angles b) differ for the Seyfert Is and the Seyfert 
2s. However, applying a KS test to the distributions puts the probability that the two sets 
of data are drawn from the same distribution as 29 percent. This is an indication that based 
on the {5, i) pairings alone (Fig. 2) one should not expect to be able to find a significant 
difference in the values of (3 for Syls and Sy2s. 

3.3 Estimation of the P(/3) distribution 

Here we describe a method whereby, in principal, we may arrive at an estimate for the 
distribution of angles (3 the jet vectors in our sample make with the galaxy disk normal. We 
denote this distribution in terms of a probability distribution P(/3), which is defined so that: 



Thus, for example if the jets axes were randomly orientated in space, and thus were inde- 
pendent of the orientation of the galaxy disk, then we would find that P(/5) = sin {3. 

If we define the jet axis in terms of the angles (5 and 6 with respect to the coordinate 
axes OXYZ, so that 



(see Fig. 1), then we note that the probability distribution P(/3) corresponds to the full 
two-dimensional probability distribution P(/9, 9) integrated over the azimuthal angle 9. In 
order to proceed further we need to make the assumption that the full two dimensional 
distribution P(/?,^) is in fact independent of 6. This assumption is a simple one and is valid 
provided that there are no selection effects which act preferentially to select in favour of jets 
with particular values of the azimuth 9. Since the angle 9 is defined in terms of coordinates 
which have a particular orientation relative to our line of sight, such a selection effect cannot 
be ruled out completely. Indeed in the Standard Model of AGN, in which the difference 
between Seyfert Is and Seyfert 2s depends on whether is greater or less than 0c, it is clear 
that this assumption is not a good one when applied to Seyfert Is or Seyfert 2s separately. 
However, when applied to a complete sample of Seyferts as a whole (which we do not have), 
the assumption may indeed be a good one. 

We begin by noting that if we were able to measure the actual value of /3 for each galaxy 
jet, then again assuming that there are no significant selection effects, we could simply plot 
a histogram of P(/9) directly from the observed sample. However we are unable to determine 
an exact value of j3 for each galaxy in our sample, because we do not know where on the 
great circle the jet vector lies - that is we do not know the value of in Equation (3.1.2). We 
can recognize, however, that if we already knew P(/5) (which we are of course attempting to 
determine) then we would be able to convert this information into a probability distribution 
for 0, that is we would be able to determine the likehhood of the jet vector lying at various 
points along the great circle. This is the basis of the method which we employ. We also note 




(3.3.1) 



= (sin (3 cos 9, sin l3 sin 9, cos 13) 



(3.3.2) 
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that of course, as discussed above, if we knew in addition the sign of and/or the sign of 
kj.ks = cos 0, then we could further restrict the position of the jet direction along the great 
circle defined by i and 6 (see Section 5). 

We start by assuming the we know P(/3), and in particular we start from the assumption 
that the jets are randomly orientated so that P(/5) = sin j3. We call this initial assumed 
distribution Pin{l3). This implies that along a great circle, the assumed local surface (i.e. 
two-dimensional) probability density of jet orientations is Pj„(/3)/ sin/?. If the jet directions 
were randomly orientated, then the probability of finding a jet vector in a particular galaxy 
at a point on the unit circle between and (f) + d(j) would be proportional to | sin0| (i0. The 
necessity for the factor | sin 0| can be seen by noting that it enables the method to reproduce 
correctly a randomly orientated jet distribution. Thus for a given Pj„(/5) the probability of 
finding the actual jet between and (p + dcj) would be [Pm(/?)/ sin/?]| sin0| dcj). 

We now proceed by taking each data point and distributing it over the unit sphere 
according to the above formula derived from the assumed Pj„(/?). Having done that we may 
then integrate over azimuth on the sphere and obtain another estimate of the (3 distribution, 
which we call PoutiP)- Clearly if and Pout are the same (or do not differ significantly), 
then we have a good estimate of the real j3 distribution. If they differ significantly, however, 
then we may, in principle, use the above method iterativelyQ 

Because we have only a small sample (even when both the Syls and the Sy2s are taken 
together) we simplify the method by binning the data. We bin the data in 5 bins equally 
spaced in cos /?, so that if the jets are randomly orientated we expect equal numbers of 
galaxies to fall into each bin. In our estimation procedure, we note that, since we redistribute 
each galaxy among the bins according to some probability distribution, the actual number of 
galaxies assigned statistically to each bin may not be a whole number. This in turn implies 
that the statistical analysis of the procedure is not straightforward. For simplicity, we adopt 
the use of -y/n error bars, where n is the (not necessarily whole) number in each bin. We 
expect this to give us a reasonable estimate of the likely error due to sampling. 

In Figure 4 we show the result obtained for Pout{P), using the input assumption that 
PiniP) corresponds to random jet orientations, or in other words, the input assumption is 
that each jet could lie anywhere randomly along the great circle defined by i and 6. As can 
be seen from the Figure, there is no evidence in the current sample that the jets orientations 
differ significantly from being uniformly distributed. We note that this result does not 
contradict the findings of Schmitt et al. (1997) who concluded, in line with our findings in 
Section 3.2, that there is a significant underabundance of AGN with (3 < 20°. Our choice of 
bin size as shown in Figure 2 means that the first bin corresponds to 0° < /3 < 37°, which 
results in such an effect being washed out in the analysis. 

3.4 The Seyfert 1 and Seyfert 2 Subsamples 

^We note briefly that the simple iterative procedure as stated here is not well posed, and that an extra 
constraint such as an assumption about the smoothness of the (3 distribution, or its closeness to being 
uniform, is likely to be required. 
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In all of the above we have not differentiated between the AGN according to whether 
they arc Scyfert Is or Seyfert 2s. However, if we are prepared to make the hypothesis that 
the Standard Model of AGN is correct, then we can make further statistical tests. We 
shall, for illustration, take the strong hypothesis that AGN are either Seyfert Is or Seyfert 
2s depending on whether is less than or is greater than (respectively) some canonical, 
and universal, value (pc- Given this assumption we may make an estimate for the P{P) 
distribution for the Seyfert Is and Seyfert 2s separately. For the Seyfert Is we assume that 
each jet can be anywhere along the great circle, with the proviso that < 0c- Similarly for 
the Seyfert 2s we make the same assumption except that now we require > 0c. 

However, using the assumption that the value of is the only difference between the 
Seyfert Is and the Seyfert 2s, it follows that the P{/3) distributions we obtain for the two 
subsamples should be the same. This means that if we vary the assumed value of 0c until the 
difference between the two P{f3) so obtained is a minimum, then we obtain an estimate of the 
value of 0c. We have carried out this calculation, using only 3 bins in cos(/3), since the Seyfert 
1 sample only contains 14 entries. We use a simple test (with two degrees of freedom) to 
see if the two distributions differ. We find that the difference between the two distributions 
is a minimum at 0c — 40°, but that the goodness of fit is poor, with the probability that the 
two distributions are different being 96 per cent (i.e approximately 2a). We show the 'best 
fit' in Figure 5. This implies that the simple hypothesis we have made in order to determine 
an estimate for 0c is not sustained. It may be that 0c is not a universal constant, for example 
it could be a function of (3, or the opening of the torus might not be cylindrically symmetric. 
Alternatively, it might be that the [3 distributions for Syl and Sy2 galaxies are significantly 
different. A possible hint in this direction comes from the observation that for the 'best 
fit' (and indeed the fits for all values of 0c) the Seyfert 1 distribution is more concentrated 
towards low values of /3 than is the case for Seyfert 2s. Since this difference between the Syl 
and Sy2 distributions is not apparent in the Pmin distributions (Section 3.2) it seems that 
it is the assumption that the Syls and Sy2s are differentiated in terms of 0, an assumption 
which is fundamental to the Standard Model, which drives the suggestion of the difference 
in the (3 distributions. 

Given the lack of systematic selection of the sample, the small sample size, and the 2a 
level of significance of the result, this result can only be called suggestive and not conclusive. 

4. Conclusions and Discussion 

We have investigated the angle j3 which the jet axis kj makes with the normal to the 
galaxy disk plane. By considering the minimum value of f3 which is allowed by the data 
for each galaxy, we are able to conclude (Section 3.2) that for both Syls and Sy2s most of 
then have angles /3 which are greater than 25° — 30". This is a straightforward conclusion, 
independent of theoretical models, and implies not only that jets perpendicular to the disk 
plane are not the norm, but also that strict alignment is avoided. 

Because the data are all derived from a projection of a three-dimensional situation onto 
the two-dimensional plane of the sky, we cannot expect to be able to reconstruct the full three- 
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dimensional picture without resorting to some statistical modelling. We have demonstrated 
how such modelling can be carried out (Section 3.3) with a minimum of assumptions. Our 
estimate of the f3 distribution for the sample as a whole (Fig. 4) shows that the distribution 
is consistent with the distribution which would result from the jet being pointed randomly 
with respect to the galaxy. 

By making the further hypothesis that 0c is a universal constant, we have been able to 
estimate a value for 0c, by choosing that value which makes the deduced j3 distributions 
for the Syls and the Sy2s most consistent with each other. We obtain the value 0c — 40°, 
but note that the two distributions so obtained for P{(5) are significantly different. Thus 
simple-minded application of the Standard Model appears to lead to a contradiction. The 
assumption that Syls and Sy2s are differentiated in terms of the angle seems to result in 
general in a Syl f3 distribution which displays significantly smaller values of f3 than the Sy2 
distribution. Using a simple test to cvahiate the significance of this result, the hypothesis 
that the Syl and the Sy2 distributions, P{P), are the same is ruled out at the 4 per cent 
level, or approximately 2a. Given the lack of systematic selection criteria of the sample, 
and the fact that similarity is ruled out only to 2a, this result is not conclusive but only 
suggestive. A complete sample, selected on an isotropic property, is required to obtain a 
conclusive result. 

Thus, we have established the feasibility of this kind of analysis, which provides ob- 
servational measures of the (3 distribution and provides some hope that the assumptions 
underlying the Standard Model for AGN might be subject to observational scrutiny. It is 
evident that a larger sample would be of value. With a large enough sample, it will also 
be necessary to take proper account of measurement errors in the angles i and S, and also 
to have a better understanding of the selection effects present in the sample. One example 
of such an effect is the probabihty that small values of are hkely to be ruled out of the 
sample, because severe foreshortening of the jet means that it is impossible to measure 5 with 
any precision. Even without a larger sample, however, further constraints on the geometry 
obtained by further observations of the sample galaxies would strengthen the conclusions. 

One such constraint would be to determine which side of the minor axis lies physically 
closer to us. If we can achieve this, we can remove the degeneracy in our knowledge of the 
sign of and in effect improve the statistics in the sample by about a factor of two since 
we can constrain the jet to lying along only a particular segment of the great circle. One 
method by which we can discover which is the nearer edge of the galaxy is by taking long slit 
spectroscopy along both the major and the minor axis. Along the minor axis we would look 
for dust effects and thus need the bluest possible waveband, where the dust has the most 
effect. By differencing the spectrum of the top minor axis from that of the bottom minor axis 
we should be able to distinguish the foreground, less-reddened side from the background, 
more reddened side. Along the major axis we would look for velocity differences, and thus 
need spectral regions with distinct galactic spectral features, which requires slightly redder 
spectral regions. Then, along with the assumption that the spiral arms are always trailing, 
the velocities along extreme ends of the major axis, or the farthest ends intercepted by a 
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spiral arm, would easily allow us to distinguish front from back. The differences in velocity 
that we will need to detect will be of order 600 km/s independent of the inclination i. By 
placing the long slit on both the major and minor axes, we determine foreground side versus 
background side with two independent methods. Both methods are more accurate for edge- 
on than for face-on galaxies, but similarly, the answer is more constraining for the edge-on 
than for the face-on case. 

A further, though more difficult, method of constraining the geometry would be to at- 
tempt to determine whether the radio jet is moving towards us as opposed to away from us 
(not to be confused with whether the radio jet is on the front side of the galaxy as opposed 
to the back side of of the galaxy) . This would then hmit the sign of cos (p and so improve our 
statistics by further limiting the arc on the great circle along which the jet can lie. Deter- 
mining whether the jet is into the sky plane or out of the sky plane is not easy. For some of 
the stronger, more highly polarized jets, it may be possible to look for depolarization along 
the length of the jet. This method is very effective for the grand radio jets, as demonstrated 
by Laing (1988) and Garrington et al. (1988) but may be problematic for the smaller jets 
found in Seyfert galaxies. As a polarized jet is observed further into its own local halo, the 
halo serves to depolarize the jet, thus showing which end of the jet is observed through a 
higher column density. 

Finally we note that although further observations, both in terms of a larger sample, 
and in terms of further observations of the current sample, will provide better constraints, 
it is clear that we have already demonstrated the power of this kind of approach. Even 
the preliminary results which we present here, based on simple observations of an initial 
sample, provide interesting implications for theoretical interpretation of Seyfert galaxies. In 
particular it may be necessary to provide a theoretical explanations (a) for why the inner 
accretion disks (as evidenced by the jets) appear to have little knowledge of the plane of the 
galaxy and (b) for why either there is little difference in the (f) values for Syls and Sy2s or 
the inner disks in Syls are more aligned with the galaxy plane than the disks in Sy2s. It will 
be of interest to see whether these results hold up with the acquisition of further and better 
observational data. 
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Figure Captions 

Figure 1. The plane of the galaxy is the X, Y plane, placed so that the apparent major 
axis of the galaxy is the X axis. The line of sight, Z , thus lies in the plane of the paper, 
and i is the angle of inclination of the galaxy. The radio jet is represented by the arrow (and 
referred to in the text as kj). The angle at which the central accretion disk axis lies relative 
to the host galaxy disk is given by 13. The angle between the line of sight and the radio axis 
is (f). In the Standard Model, active galaxies with (f) less than some canonical value of about 
30° are observed as Scyfcrt 1 galaxies, while active galaxies with (f) greater than that value 
are observed to be Seyfert 2 galaxies. The position angle between the radio axis and the 
apparent major axis of the galaxy is given by 5- Note that the two observable angles are i 
and 5. The azimuthal angle of the radio jet is given by 9. 

Figure 2. The difference in position angle on the sky (dpa) between the jet and the apparent 
major axis of the galaxy (also referred to in the text as 5) is plotted against the inclination 
angle i of the galaxy plane (with i — 90° corresponding to an edge-on galaxy) . The Seyfert 
Is are shown as filled squares and the Seyfert 2s as crosses. 

Figure 3. The cumulative distributions of I3min are shown for both the Syls and the 
Sy2s. The distributions are scaled so that the total number of each subsample appear as a 
percentage. The Seyfert Is are shown as filled squares and the Seyfert 2s as crosses. 
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Figure 4. The estimated P{P) distribution obtained for the whole sample (Syls and Sy2s). 
The bins are chosen to have equal sizes in 1 — cos{(3) so that a uniform distribution on the 
sky would give equal numbers in each bin (as shown by the horizontal line). The error bars 
in each bin have been simply estimated as ^/n where n is the number assigned to the bin. 
Figure 5. The estimates of the P{(3) distributions for (a) the Syls and (b) the Sy2s which 
result from the assumption that the Syls and Sy2s are differentiated only in terms of whether 
(f) is less than, or greater than, the 'best fit' value of 0c — 40°. As for Figure 4 the bin sizes 
are equal in 1 — cos{(3). 
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